Perturbation of multiparameter 
non-self-adjoint boundary eigenvalue 
problems for operator matrices 

O.N. Kirillov 

Abstract. We consider two-point non- self- adjoint boundary eigenvalue prob- 
lems for linear matrix differential operators. The coefficient matrices in the 
differential expressions and the matrix boundary conditions are assumed to 
depend analytically on the complex spectral parameter A and on the vector of 
real physical parameters p. We study perturbations of semi-simple multiple 
eigenvalues as well as perturbations of non-derogatory eigenvalues under small 
variations of p. Explicit formulae describing the bifurcation of the eigenvalues 
are derived. Application to the problem of excitation of unstable modes in ro- 
tating continua such as spherically symmetric MHD a 2 -dynamo and circular 
string demonstrates the efficiency and applicability of the theory. 
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1. Introduction 

Non-self-adjoint boundary eigenvalue problems for matrix differential operators 
describe distributed non-conservative systems with the coupled modes and appear 
in structural mechanics, fluid dynamics, magnetohydrodynamics, to name a few. 

Practical needs for optimization and rational experiment planning in modern 
applications allow both the differential expression and the boundary conditions to 
depend analytically on the spectral parameter and smoothly on several physical 
parameters (which can be scalar or distributed). According to the ideas going 
back to von Neumann and Wigner [T], in the multiparameter operator families, 
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eigenvalues with various algebraic and geometric multiplicities can be generic |12j . 
In some applications additional symmetries yield the existence of spectral meshes 
[33| in the plane 'eigenvalue versus parameter' containing infinite number of nodes 
with the multiple eigenvalues [TU1 E3 [33 ISH] • As it has been pointed out already by 
Rellich [5] sensitivity analysis of multiple eigenvalues is complicated by their non- 
differentiability as functions of several parameters. Singularities corresponding to 
the multiple eigenvalues [12] are related to such important effects as destabilization 
paradox in near-Hamiltonian and near-reversible systems HH [25l [26l [271 EH] > 
geometric phase [31], reversals of the orientation of the magnetic field in MHD 
dynamo models [32] . emission of sound by rotating continua interacting with the 
friction pads [38] and other phenomena [23] . 

An increasing number of multiparameter non-self-adjoint boundary eigen- 
value problems and the need for simple constructive estimates of critical parame- 
ters and eigenvalues as well as for verification of numerical codes, require develop- 
ment of applicable methods, allowing one to track relatively easily and conveniently 
the changes in simple and multiple eigenvalues and the corresponding eigenvec- 
tors due to variation of the differential expression and especially due to transition 
from one type of boundary conditions to another one without discretization of the 
original distributed problem, see e.g. [TT1 [15| [T8 ] [20 ] [36 ] [37 ] [38] . 

A systematical study of bifurcation of eigenvalues of a non-self-adjoint linear 
operator Lq due to perturbation Lq + eLi, where e is a small parameter, dates back 
to 1950s. Apparently, Krein was the first who derived a formula for the splitting of a 
double eigenvalue with the Jordan block at the Hamiltonian 1 : 1 resonance, which 
was expressed through the generalized eigenvectors of the double eigenvalue [3] . In 
1960 Vishik and Lyusternik and in 1965 Lidskii created a perturbation theory for 
nonsymmetric matrices and non-self-adjoint differential operators allowing one to 
find the perturbation coefficients of eigenvalues and eigenfunctions in an explicit 
form by means of the eigenelements of the unperturbed operator [H [S] . Classical 
monographs by Rellich b 9j , Kato [7J , and Baumgartel [T3] , mostly focusing on the 
self-adjoint case, contain a detailed treatment of eigenvalue problems linearly or 
quadratically dependent on the spectral parameter. 

Recently Kirillov and Seyranian proposed a perturbation theory of multiple 
eigenvalues and eigenvectors for two-point non-self-adjoint boundary eigenvalue 
problems with scalar differential expression and boundary conditions, which de- 
pend analytically on the spectral parameter and smoothly on a vector of phys- 
ical parameters, and applied it to the sensitivity analysis of distributed non- 
conservative problems prone to dissipation-induced instabilities [HI [T7J HH [211 
[20l l24l |26"1 l27l [38] . An extension to the case of intermediate boundary conditions 
with an application to the problem of the onset of friction-induced oscillations in 
the moving beam was considered in |34j . In |33j this approach was applied to the 
study of MHD a 2 -dynamo model with idealistic boundary conditions. 

In the following we develop this theory further and consider boundary eigen- 
value problems for linear non-self-adjoint m-th order N x N matrix differential 
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operators on the interval [0, 1] 3 x. The coefficient matrices in the differential ex- 
pression and the matrix boundary conditions are assumed to depend analytically 
on the spectral parameter A and smoothly on a vector of real physical parameters p. 
The matrix formulation of the boundary conditions is chosen for the convenience of 
its implementation in computer algebra systems for an automatic derivation of the 
adjoint eigenvalue problem and perturbed eigenvalues and eigenvectors, which is 
especially helpful when the order of the derivatives in the differential expression is 
high. Based on the eigenelements of the unperturbed problem explicit formulae are 
derived describing bifurcation of the semi-simple multiple eigenvalues (diabolical 
points) as well as non-derogatory eigenvalues (branch points, exceptional points) 
under small variation of the parameters in the differential expression and in the 
boundary conditions. Finally, the general technique is applied to the investigation 
of the onset of oscillatory instability in rotating continua. 



2. A non-self-adjoint boundary eigenvalue problem for a matrix 
differential operator 

Following [8] [19l [22l [24l [26], [27] we consider the boundary eigenvalue problem 

L(A,p)u = 0, U fe (A,p)u = 0, k=l,...,m, (1) 

where u(x) GC W « C( m '[0, 1]. The differential expression Lu of the operator is 

m 

Lu = ^lj(a;)d™-Ju, l 3 {x) G C NxN <g> C (m - j) [0, 1], detpo(ar)] + 0, (2) 

3=0 

and the boundary forms U^u are 

m — 1 m—1 

U k u = J2 A fe ,u(% = 0) + £ B kj u^(x = 1), A fe „ B kj G C™. (3) 

3=0 3=0 

Introducing the block matrix it := [21, 93] G C mNx2mN and the vector 

u T := (u r (0), u« T (0), . . . , u^-^CO), u T (l), u^ T (l), u^ T (l)) G C 2 "^ 

(4) 

the boundary conditions can be compactly rewritten as [261 127j 

Hu = [21, «8]u = 0, (5) 

where 21 = (A kj )\ x=Q G C mNxmN and OS = (B kj )\ x=1 G C mNxmN . It is assumed 
that the matrices lj, 21, and S are analytic functions of the complex spectral 
parameter A and smooth functions of the real vector of physical parameters p € M. n . 
For some fixed vector p = p the eigenvalue Ao, to which the eigenvector Uo 
corresponds, is a root of the characteristic equation obtained after substitution of 
the general solution to equation Lu = into the boundary conditions ([5]) [5]. 
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Let us introduce a scalar product [8] 



<u,v>:= / v*uc?x, (6) 
Jo 

where the asterisk denotes complex-conjugate transpose (v* := v T ). Taking the 
scalar product of Lu and a vector-function v and integrating it by parts yields the 
Lagrange formula for the case of operator matrices (cf. [H [22] [26] [27] ) 

Sl(u, v) :=< Lu, v > - < u,L f v >= 0*£u, (7) 

with the adjoint differential expression [5] [55] 



Ltv:=^(-1)™-^™-" (l*v) 



(8) 



the vector D 



» T := (/(O),v( I ) ^ (0),...,v^^(0),/(l),vf(l),..,v^ I ) ^ (l))GC M 

(9) 

and the block matrix C := (ly) 



C 



-£(0) 




£(1) ) ' £(x) 



/ loo 
lio 



loi 

111 



V i 



lm-20 lm-21 
?n- 10 



lom-2 lom-1 \ 
llm-2 





o / 



(10) 



where the matrices Uj are 



m — l—j 

(-l) k M% d^lm-l-j-k, 

k—i 

' ' i + j < m - 1 n fc > i > 



M 




(11) 



Extend the original matrix it (cf. ([5])) to a square matrix which is made 
non-degenerate in a neighborhood of the point p = Po and the eigenvalue A = Ao 
by an appropriate choice of the auxiliary matrices 2t(A,p) and *8(A,p) 

21 S 



it = [21, <B] ^ U := 



21 <B 



i2mNx2mN 



, il:= [21, <8], det(U) ^ 0. (12) 



Similar transformation for the adjoint boundary conditions 23 = [£, D]v = 
yields 



2J := [£,£] ^ V := 



^2rnN x2mN 



, 23 := [£,£>], det(V) ^ 0. (13) 
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Then, the form fi(u, v) in ([7]) can be represented as [8] 

0(u,v) = (9Jt))*Ik - (2Jt))*Uu, 
so that without loss in generality we can assume [551 HI] 

C = 27*H - 2TIt. 
Differentiating the equation (|T5|) we find 



k=0 



(14) 
(15) 

(16) 



Hence, we obtain the formula for calculation of the matrix 2J of the adjoint bound- 
ary conditions and the auxiliary matrix 2J 



-aj 

53 



£(1) 



21 <B 



(17) 



which exactly reproduces and extends the corresponding result of [26[ [2?1 



3. Perturbation of eigenvalues 

Assume that in the neighborhood of the point p = Po the spectrum of the bound- 
ary eigenvalue problem |T|) is discrete. Denote L =L(Ao,Po) and ilo=il(Ao, Po). 
Let us consider a smooth perturbation of parameters in the form p = p(e) where 
p(0) = po and e is a small real number. Then, as in the case of analytic matrix func- 
tions [55], the Taylor decomposition of the differential operator matrix L(A, p(e)) 
and the matrix of the boundary conditions it(A, p(e)) are [HO [2H HSj [57] 

L(A,p(e)) = f) ^e-U U(A,e) = £] ^-^L e ^ rs , (18) 

r,s— r,s— 

with 

Loo = Lo, L r o = d^L, 

i^oo = ilo; ^rO = c^il, 

where dot denotes differentiation with respect to e at e = and partial derivatives 
are evaluated at p = po, A = Ao. Our aim is to derive explicit expressions for 
the leading terms in the expansions for multiple-semisimple and non-derogatory 
eigenvalues and for the corresponding eigenvectors. 



(19) 
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3.1. Semi-simple eigenvalue 

Let at the point p = po the spectrum contain a semi-simple /i-fold eigenvalue 
Ao with n linearly- independent eigenvectors Uo(x), Ui(cc), . . ., u Al _i(x). Then, the 
perturbed eigenvalue A(e) and the eigenvector u(e) are represented as Taylor series 
me [! [MUSI US E3] 

A = A +eAi +e 2 A 2 + . . . , u = b + ebi + e 2 b 2 + . . . . (20) 

Substituting expansions (fT8"|) and ([217)1 into ([T]) and collecting the terms with 
the same powers of e we derive the boundary value problems 

L b - 0, il b - 0, (21) 

L bi + (AiLi + L i)b = 0, ilobi + (Aiilio + H i)b = 0, (22) 

Scalar product of (f2"2")l with the eigenvectors v,-, j = 0, 1, . . . , fi — 1 of the adjoint 
boundary eigenvalue problem 

L^v = 0, 2J = (23) 

yields /i equations 

< Lobi,Vj >= - < L 01 b ,v, > - Ai < Liob ,v,- > . (24) 

With the use of the Lagrange formula ([7]) , (TT5]) and the boundary conditions (12"2")) 
the left hand side of (|24|) takes the form 

< L bi, v, >= (» o,)*(il ibo + Aiil 10 b ). (25) 
Together ([24]) and ([25]) result in the equations 

Ai (< Lxobo.Vj > + (» t>i)*H l0 bo) = - < L 01 b ,v, > - (W t>i)*iIoil»o. (26) 

Assuming in the equations (|26p the vector bo(x) as a linear combination 

bo (a;) = c u (a;) + ciUi(x) + . . . + c /1 _iu /1 _i(x), (27) 

and taking into account that 

b = c Uo + ciui + . . . + fy-iu^-i, (28) 

where c T = (cq, ci, . . . , c^-i), we arrive at the matrix eigenvalue problem (cf. [?]) 

Fc = AiGc. (29) 

The entries of the /i x /i matrices F and G are defined by the expressions 

Fij =< L Q1 u i ,v i > + D*2Jol"toi%i Gij =< L 10 Uj,Vi > + o*^S^io%. (30) 

Therefore, in the first approximation the splitting of the semi-simple eigenvalue 
due to variation of parameters p(e) is A = Ao + eAi + o(e), where the coefficients 
Ai are generically \x distinct roots of the /z-th order polynomial 

det(F + AiG) = 0. (31) 
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For n = 1 the formulas ([30]) and (|3Tj) describe perturbation of a simple eigenvalue 

A = Ao — £ — 1- o(e). (32) 

< Li u ,v > + 05*05 iiiouo 

The formulas ([30]) , ([3~T|) , and (|32p generalize the corresponding results of the works 
[24| to the case of the multiparameter non-self-adjoint boundary eigenvalue prob- 
lems for operator matrices. 

3.2. Non-derogatory eigenvalue 

Let at the point p = Po the spectrum contain a /z-fold eigenvalue Ao with the 
Keldysh chain of length /i, consisting of the eigenvector Uo(x) and the associated 
vectors Ui(x), . . .,u^_x(ai). The vectors of the Keldysh chain solve the following 
boundary value problems [5J [22] 

L u = 0, H Uo = 0, (33) 
1 I j 1 

L Uj = - 22 — ^A Lu i-r. i^oUj = - X/ — 5 A^Uj- r . (34) 
r— 1 r— 1 

Consider vector-functions Vo(a;), vi(x), . . ., v^_i(a;). Let us take scalar prod- 
uct of the differential equation ([3"3")l and the vector- function v M _i(a;). For each 
j = 1, . . . , fi — 2 we take the scalar product of the equation (|34[) and the vector- 
function v^-i-j (.t). Summation of the results yields the expression 



EE h < 9 I Lu i-> v M-i-i >= (35) 

Applying the Lagrange identity (J7J), (fT5|) and taking into account relation (flTjj) . we 
transform (13511 to the form 



E < u c-H^^ Lfv i-' > 

j=Q r=0 

/x— 1 fj,— X—k 

E E 

fc=0 3=0 



r=0 



-j^-x-j-* = 0. (36) 



Equation (|36p is satisfied in case when the vector- functions vo(x), Vx(x), . . ., 
v fl —i(x) originate the Keldysh chain of the adjoint boundary value problem, cor- 
responding to the //-fold eigenvalue Ao [26l [27] 



lJv - 0, 2J D - 0, (37) 

j J 

L ^ = - E ^ L Vr, 5J 0, = - E i^^-r- (38) 



r— 1 r— 1 
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Taking the scalar product of equation (|34[) and the vector Vo and employing 
the expressions ([7]), (TT5)) we arrive at the orthogonality conditions 

$3 - [< ^Lu 3 - r , vo > + t^QJ^Hu^J =0, j = 1, . . . , M - 1. (39) 

r=l 

Substituting into equations (fl]) the Newton-Puiseux series for the perturbed 
eigenvalue A(e) and eigenvector u(e) [7J H21 [57] 

A = A + Aie 1 /' 1 + . . . , u = w + w^ 1 /" + . . . , (40) 

where w = u , taking into account expansions (|18[) and (j4"0)) and collecting terms 
with the same powers of e, yields /i — 1 boundary value problems serving for 
determining the functions wy, r — 1,2, ... ,fj,— 1 



L w r = I ^ -^L ff0 J! A Ql ...A Q ] Wj , (41) 

3=0 \ct=1 ' \a\„=r—j ) 

r-lr-j / \ 1 

U tDr.--XX X, ^•■■ A «J^0tl}, (42) 

where |a| CT = ai+ . . . +a cr and a\, . . ., are positive integers. The vector- 

function w M (x) is a solution of the following boundary value problem 

M-i / A 
Low M = -LqiWq - X X — L ct0 X A ai ...A a I w.,, (43) 
i=o \ct=i ' H^/i-j / 

p-i/j-j / \ , 

Hotn M = -iloiWo - A Q1 . . . A Q I —ii^Wj. (44) 

3=0 cr=1 \|os|<r=/*— J / 

Comparing equations (f4"3")l and (|4~4"1) with the expressions (j3~4"|) we find the first fi— 1 
functions w r in the expansions (|40[) 

r 

w r = X Uj X A Ql . . . \ aj . (45) 

3=1 |a| 3 =r 

With the vectors (|4"S")| we transform the equations (|4*3"|) and (|4"4")l into 

L w M = -Loiu - A^^— ^Lu M _ r + ^ L Uj ^ A ai ...Aa 3 , (46) 

r-=l ' 3=1 |«b=A« 

itoW^ = -UiUo- j^-j'^HUp-r + J^Hott,- X A ai---^a r (47) 



3=1 |a|j=M 
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Applying the expression following from the Lagrange formula 



< L w„, v > 



M 1 _ 

o*Qj*iioiu + x^22 -f&sajjUrou,,- 



(48) 



and taking into account the equations for the adjoint Keldysh chain (|37[) and 
yields the coefficient Ai in (|40"]) . Hence, the splitting of the /i-fold non-derogatory 
eigenvalue Ao due to perturbation of the parameters p = p(e) is described by the 
following expression, generalizing the results of the works [HI [24J [26j [27] 



A = An + a 



< L iUo,v > + 



+ o(e" ). 



For /i = 1 equation (|49|) is reduced to the equation (|32[) for a simple eigenvalue 



(49) 



4. Example 1: A rotating circular string 

Consider a circular string of displacement W(ip,r), radius r, and mass per unit 
length p that rotates with the speed 7 and passes at <p — through a massless 
eyelet generating a constant frictional follower force F on the string, as shown in 
Fig. [1] The circumferential tension P in the string is assumed to be constant; the 
stiffness of the spring supporting the eyelet is K and the damping coefficient of the 
viscous damper is D. Introducing the non-dimensional variables and parameters 




and assuming w(<p,t) = u(ip)exp(Xt) we arrive at the non-self-adjoint boundary 
eigenvalue problem for a scalar (N = 1) differential operator [T5] 

Lu = \ 2 u + 2n\u' - (l-n 2 )u" = 0, (51) 

U (0)-u(2tt) = 0, u'(0)-u'(2n) = T -^u(0) + T -^u\0), (52) 

where prime denotes differentiation with respect to ip. Parameters J7, d, fc, and fi 
express the speed of rotation, and damping, stiffness, and friction coefficients. 

For the unconstrained rotating string with d = 0, k = 0, and /x = the eigen- 
functions v and u of the adjoint problems, corresponding to purely imaginary eigen- 
value A and A, coincide. With u = C\ exp (<p\/(l — O)) + C2 exp (—ip\/(l + fi)) 
assumed as a solution of ([5T|) . the characteristic equation follows from (T52"|) 
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ImA. 





Figure 1. A rotating circular string and its spectral mesh (only 
30 modes are shown). 



The eigenvalues A^ = in(l ± fi) with the eigenfunctions = cos(mp) =p 
i sin(ncp) , n € Z, found from equation (|53[) form the spectral mesh in the plane 
(fi,ImA), Fig.Q] The lines A^ = in(l + eSl) and A^ = im(l + 5SY), where e, 5 = ±1, 
intersect each other at the node (0^ n , A^„) with 



n 



eS 



X 



eS 



inm{5 - 



(54) 



(55) 



mo — ne too — ne 

where the double eigenvalue X^ n has two linearly independent eigenfunctions 

u 6 n = cos(rup) — ei sin(nip), u s m = cos(rrnp) — Si sin(my). 

Using the perturbation formulas for semi-simple eigenvalues ()30|) and (|3ip 
with the eigenelements (f54| and (|55p we find an asymptotic expression for the 
eigenvalues originated after the splitting of the double eigenvalues at the nodes 
of the spectral mesh due to interaction of the rotating string with the external 
loading system 



nra 



.en 



8m . _ n 

— An + % 



87mm 



(dX^t + k) + 



8- 



-M±\/c, 



(56) 



with An = Q - nfL and 



en — 8m ^ m — n , 

■Afi + i- (dX 



(dX E n A m + k - ien f i){dX £ n d m + k 



8nmn 

eS 



6 
tun 



+ k) + 
iSmfj,) 



8tt 



1Q-K 2 nm 



(57) 



nrn 



Due to action of gyroscopic forces and an external spring double eigenvalues 
split in the subcritical region |f2| < 1 (e < 0, <5 > and to > n > 0) as 



x = x±+ l T ^An 



.n + m. , . 

+ 1 K±l\ 



k 2 



/to — 



to + n 



16n 2 nm V 8Trmn 



Afi , (58) 
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Figure 2 . Deformation of the spectral mesh of the rotating string 
interacting with the external spring with k = 2 (left); approxima- 
tion to the corresponding tongues of oscillatory instability (right). 



while in the supercritical region |0| > 1 (e < 0, 5 > and m > 0, n < 0) 



.. s .m + \n\ . ^ ,\n\—m, , 



k 2 



87r|n|m 



167T 2 |7l|m 



-Aft 



87rm|? 



(59) 

Therefore, for |f2| < 1 the spectral mesh collapses into separated curves demon- 
strating avoided crossings; for |f2| > 1 the eigenvalue branches overlap forming the 
bubbles of instability with eigenvalues having positive real parts, see Fig. [2] From 
(jB"9"]) a linear approximation follows to the boundary of the domains of supercritical 
flutter instability in the plane (f2, k) (gray resonance tongues in Fig. [2]) 

47r|n|m(|n| — m) 



— . o - f-H — • (6°) 

At the nodes with fl^ m = the external damper creates a circle of complex 
eigenvalues being a latent source of subcritical flutter instability responsible for 
the emission of sound in the squealing brake and the singing wine glass [38 

2 d 2 
ReA + — ) + n 2 il 2 = — — T , ImA = n, 



n 2 n 2 - (ImA — nY 



16tH 



ReA 



d_ 

"4tt' 



16tt 2 ' 

as shown in Fig. [3l Non-conservative perturbation yields eigenvalues with 

1 



2n 



ImA = n ± — V 2ir 2 n 2 Q 2 ± TrnQy '4ir 2 n 2 n 2 



ReA = ±- 



1 

"2^ 



-2ir 2 n 2 n 2 ± TtnVL^A-K 2 n 2 n 2 + /i 2 , 



(61) 
(62) 

(63) 
(64) 
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Figure 3. Deformation of the spectral mesh of the string near 
the nodes (0,3), (0,2), and (0,1). 



so that both the real and imaginary parts of the eigenvalue branches show a 
degenerate crossing, touching at the node (0, n), Fig. [3] Deformation patterns of the 
spectral mesh and first-order approximations of the instability tongues obtained 
by the perturbation theory and shown in Fig. [2] and Fig. [3] are in a good qualitative 
and quantitative agreement with the results of numerical calculations of [15] . 



5. Example 2: MHD a 2 -dynamo 

Consider a non-self-adjoint boundary eigenvalue problem (N — 2) appearing in 
the theory of MHD a 2 -dynamo [H [Ml [31 [37] 



Lu := hd^u + hd x u + l 2 u = 0, Uu : 
with the matrices of the differential expression 



[21, <8]u = 0, 



In 



1 

-a(x) 1 



li = <9 X 1 , h = 



a(x) 



a{x) 
'O+i) _ A 



and of the boundary conditions 

/ 1 \ 
10 



21 = 




\ / 



<B = 



31 






- 1 





\ 



Q 

10 0/ 



(65) 
(66) 

(67) 
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where it is assumed that a(x) = ao + 7Aa(a;) with Aa(x)dx — 0. For the 
fixed Aa(x) the differential expression depends on the parameters ao and 7, while 
(3 interpolates between the idealistic ((3 — 0) and physically realistic ((3 = 1) 
boundary conditions [TH [32 [33] [36] . 

The matrix 27 of the boundary conditions and auxiliary matrix 23 for the 
adjoint differential expression L^v = lg<9 2 v — I*^xV + l£v follow from the formula 
(fT7|) where the 4x4 matrices 21, and 23 are chosen as 

a-(s;), s-(s;). ( «, 

In our subsequent consideration we assume that I = and interpret [3 and 7 
as perturbing parameters. It is known [33) that for (3 — and 7 = the spectrum 
of the unperturbed eigenvalue problem (|65p forms the spectral mesh in the plane 
(ao, A), as shown in Fig. [4] The eigenelcmcnts of the spectral mesh are 

A^j = — (7m) 2 + eccQTrn, Af„ = — (nm) 2 + Saoitm, e, 6 = ±1, (69) 
< = I * ] sin(nTrr), uf n = [ 1 ) sin(mTrr). (70) 



eirn J m \ Snm 

The branches (|69p intersect and originate a double semi-simple eigenvalue with 
two linearly independent eigenvectors (|70p at the node (c^Ag), where ;33j 

Ag = e57r 2 nm, = enn + Snm. (71) 

Taking into account that the components of the eigenfunctions of the adjoint 
problems are related as V2 = u\ and v\ = 112, we find from equations (|30|) and 
(|31|) the asymptotic formula for the perturbed eigenvalues, originating after the 
splitting of the double semi-simple eigenvalues at the nodes of the spectral mesh 

A = Aq -eSTr 2 mn(3+ ^(6m + en)Aa (72) 



± 

where 



^y / (((5m-en)Aa ) 2 + 4rnn(e7Aa-(-l)"+ m 7rn / 9)((57Aa-(-l) ,l + m 7rm/3), 

A«o := aQ — ao, Aa := / Aa(x) cos((en — 8m)irx)dx. (73) 

Jo 

When 7 = and Aao — 0, one of the two simple eigenvalues (|73J) remains unshifted 
in all orders of the perturbation theory with respect to the parameter (3: A = Ag. 
The sign of the first-order increment to another eigenvalue A = Ag — 2Ag/3 depends 
on the sign of Ag , which is directly determined by the Krein signature of the modes 
involved in the crossing [33] . This is in the qualitative and quantitative agreement 
with the results of numerical calculations of [35] shown in Fig. 0] 

Under variation of the parameter [3 in the boundary conditions the eigenval- 
ues remain real. Additional parameter 7 is required for the creation of complex 
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Figure 4. 1 = 0: Deformation of the spectral mesh for 7 = and 
(3 G (0, 1); approximation of the resonance tongues for Aa(x) = 
cos(47ra;) and Aq < (white, light gray) or Aq > (dark gray). 



eigenvalues, which happens when the radicand in (f?2")) becomes negative 

((en- 5m)Aa ) 2 + ran ((e + 5)jAa - (-l) n+ro (n + m)(3n) 2 

< mn ((s - S)jAa - (-l)"-" l (n - to)/?tt) 2 (74) 

Inequality (|74|) defines the inner part of a cone in the space (ao,(3,j). The 
part of the cone corresponding to ReA > (oscillatory dynamo) is selected by the 
condition 

Aq - eSTT 2 mn/3 + ^(Sm + en)Aa > 0. (75) 

The conical zones develop according to the resonance selection rules discov- 
ered in [33]. For example, if Aa(x) = cos(27rfca:), then 

Aa = ( cos(2irkx) cos((en — 8m)irx)dx — { y\/_! , ^ OT r (76) 



0, 2k^±(en-5m) 

There exist 2k — 1 cones for Aq < (e5 < 0) and infinitely many cones for Aq > 
(e5 > 0). Due to different inclinations of the cones, only cross-sections of 2k — 1 
cones with ReA < survive in the plane (3 — 0. They are situated symmetrically 
with respect to the 7-axis 



(Acq) 2 < £ 



k,k + l,...,2k-l. (77) 



For k = 2 three resonant tongues 4ckg < j 2 and 16 (cto ± 2w) 2 < 3j 2 are shown 
white in Fig. |H When (3^0 the tongues (|77| . corresponding to Aq < 0, are 
deformed into the domains in the plane (00,7) bounded by the hyperbolic curves 

4fc 2 (a ± a v ) 2 + ((n - m) 7 /2 ± 7r(n 2 + m 2 )[3) 2 < k 2 (7 ± 2ir(n - m)(3) 2 . (78) 

The deformed resonant tongues are located at a distance 70 = 7(Aao = 0) from 
the ap-axis, which for the tongues with afj > is greater (70 = 2-nnff) than for 
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those with ckq < (70 = 27r(2fc — n)/3). In case when k = 2 the approximation to 
the deformed principal resonant tongues 

7 2 -4a 2 , > 16tt 2 /3 2 , 16(a ± 2tt) 2 + (7 ± IOtt/3) 2 < 4( 7 ± 4tt/3) 2 (79) 

is shown light gray in Fig. 2] 

Cross-sections by the plane (3 ^ of the cones, corresponding to Ag > 0, have 
the form of the ellipses, shown dark gray in Fig. 0] 

4fc 2 (a ± a^f + n(2k +n) (7 ± 2(n + k)nf3) 2 < n{2k + n)4k 2 iT 2 (3 2 , (80) 

where n = 1,2,.... The eigenvalues inside the ellipses have positive real parts, 
which corresponds to the excitation of the oscillatory dynamo regime. The ellipses 
belong to a corridor bounded by the lines 27 = k/3(ao ±47r). Hence, the amplitude 
7 of the resonant perturbation of the a-profile ~/Aa(x) is limited both from below 
and from above in agreement with the numerical findings of [32] . 

When j3 —* 0, the ellipses shrink to the diabolical points (ag, 0) in the plane 
(ao, 7). The reason for this effect is the inclination of the cone (|50|) . The rightmost 
picture of Fig. 0] shows a cross-section of the inclined cone (f5U|) (dark gray) by the 
plane ag = together with that of the cone (|75|) (light gray) . Obviously, variation 
of 7 (or, equivalently, of a-profilc) excites the complex eigenvalues only near the 
diabolical points with Ag < in accordance with [33] , while the variation of (3 does 
not produce the complex eigenvalues at all [36]. Nevertheless, the variation of (3 
together with 7 yields the complex eigenvalues near the nodes of the spectral mesh 
with both Ag < and Ag > 0. We note that the analytical results are confirmed 
both qualitatively and quantitatively by the Galerkin-based numerical simulations. 

Conclusion 

A multiparameter perturbation theory for non-self-adjoint boundary eigenvalue 
problems for matrix differential operators is developed in the form convenient for 
implementation in the computer algebra systems for an automatic calculation of 
the adjoint boundary conditions and coefficients in the perturbation series for 
eigenvalues and eigenvectors. The approach is aimed to the applications requiring 
frequent switches from one set of boundary conditions to another. Two studies of 
the onset of instability in rotating continua under symmetry-breaking perturba- 
tions, demonstrate the efficiency of the proposed approach. 
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